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Abstract 

We study the production characteristics of 557 pairs of charmed hadrons pro- 
duced in tt~ — Cu interactions at 230 GeV/c using a momentum estimator for 
charmed hadrons with missing decay products. We find, the mean value of the trans- 
verse momentum squared of the charmed pairs is <p|. sum >= (1.98 ± 0.11 ± 0.09) 
GeV 2 /c 2 , the mean rapidity difference is < \yaiff\ >= 0.54 ± 0.02 ± 0.24, and the 
mean effective mass is <M eff >= (4.45 ± 0.03 ± 0.13) GeV/c 2 . Comparing these 
results with the next-to-leading order QCD predictions we find an agreement for 
the ydi f / and M e j f , whilst the measured mean value of p\ sum is significantly larger 
than the predicted value. 



(Submitted for publication in Physics Letters B) 



1 Introduction 



Experimental results on charm hadronic production are usually compared with the 
next-to-leading ordeiQ QCD calculations Let us recall that the ACCMOR0 and 
the E769|3J collaborations have shown that the charm production cross sections, as well 
as the xp and spectra of charmed particles, do not depend strongly on the nature of 
the charmed particles^ under consideration (D°, D + , D* + , and A+). This indicates 
that the main features of charm production are determined at the parton level with a 
relatively small dependence on the light quarks forming a charmed hadron. On the other 
hand, the above characteristics are not very sensitive to the difference between the leading 
orderf] and 0(a|) QCD calculations. The charm pair correlations are more sensitive to 
this difference. 0(a 2 s ) QCD predicts that the heavy quark pair is produced exactly in a 
back-to-back configuration, corresponding to an azimuthal angle difference of A0=18O° 
and p\ sum =0. The third order terms cause the broadening of these distributions. Thus, 
the distributions of the charmed pairs test the importance of the next-to-leading order 
terms. On the other hand, disagreement with the perturbative results may show the 
importance of various nonperturbative effects, which are likely to play an important role 
at the energy scale set by the charmed quark mass. 

There is a relative abundance of data on the azimuthal angle difference distribution 
A0 of a charmed pair in hadroproduction as this angle can be determined from the 
directions of the charmed particles without the momentum measurements, which demand 
the detection of all decay products. The results of the WA75|], E653§, ACCMORfTJ, 
and to some extent, of the WA92|| collaborations show the back-to-back enhancement 
of the A<ft distribution to be much weaker than the 0(a^) predictions. There are less 
experimental results available for distributions of other kinematical variables of a charmed 
pair because they require a knowledge of the momenta of the charmed particles, while 
the fully reconstructed decays (no missing decay products) are much less numerous. The 
E653|| collaboration attempted to overcome this difficulty by using a momentum 
estimator. In this letter we perform a similar analysis of the ACCMOR data with the 
aim of determining the following characteristics of the pair of charmed particles: 
Prsum ~ sc L uare °f the vector sum of transverse momenta; 
Udiff - the rapidity difference between particles; 
M e ff - effective mass of the pair. 

Since the ACCMOR experiment, data processing and selection of double charm events 
have already been described (see ref. J7J and references therein), we recall here only the 
essential features of the experiment and the data analysis. Next, we discuss the momentum 
estimator and its errors. Finally, the results are presented and compared to the 0(af) 
QCD predictions. 



1 Hereafter called O(al). 

2 Throughout the paper a particle symbol stands for particle and antiparticle, i.e., any reference to a 
specific state implies the charge-conjugate state as well. Thus, e.g. D + D° stands for D~D° as well. 
3 Hereafter called 0(a 2 s ). 
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2 Experiment, data analysis and acceptance correc- 
tions 



The second phase of the NA32 experiment was performed at the CERN-SPS using 
a negative 230 GeV/c beam (96% pions and 4% kaons) and a 2.5mm Cu target. Charm 
decays were reconstructed with an improved silicon vertex detector and a large acceptance 
spectrometer. The latter consisted of two magnets, 48 planes of drift chambers and three 
multicellular Cherenkov counters allowing 7T, K, p identification in a wide momentum 
range. The vertex detector consisted of a beam telescope (seven microstrip planes) as well 
as a vertex telescope with two charge-coupled devices (CCDs) and eight microstrip planes. 
The high precision vertex detector allowed the clean reconstruction of charm decays with 
very few background events and a purely topological charm search, which was restricted 
neither to a limited number of decay modes nor to any mass window. 

Events with the primary vertex inside the target and at least two secondary tracks 
not originating from the vertex were selected for further analysis. These tracks were then 
used as a seed to search for one or more secondary vertices. Events with two secondary 
vertices were next selected. Since most of them involved unseen neutral decay products, 
one could not demand the effective mass of charged decay products to be compatible with 
the mass of a charmed particle. Secondary interactions were excluded by demanding a 
separation of the secondary vertex of at least two standard deviations from the target 
edge as well as from the edges of both CCDs. The decays of strange particles were 
rejected by demanding that the effective mass m V i S of the charged decay products for 
7T + 7r - vertices had to be larger than the kaon mass. Similarly, all three-prong vertices, 
which could be interpreted as an accidental overlap of a K s — > 7r + 7r~ (or A 
decay and a track, were rejected by checking the effective mass of the 7t + tt~ (or pn~) 
combinations. Additionally, the visible transverse momentum p^ s , with respect to the 
direction of the parent charmed particle P, was demanded not to exceed the maximum 
transverse momentum kinematically allowed for any particle in the decay channel under 
consideration. Finally, decays attributed to short-lived D°, or A+ had to occur before 
the second CCD (20 mm from the target). 

In ref. |l0l such selected secondary vertices were assigned to various decay modes of 
D°, D + , Df and A+. This was done with the help of m V i S and of the neutral mass defined 
as 

,2 _o m _./3 



m =m P + m vis - 2m P ^m 2 vis + (p v T ls ) 2 . (1) 

For the purpose of this letter the exact assignment is not fully needed since we only want 
to identify the charmed particle P, as a knowledge of the decay channel is irrelevant here. 
To use the momentum estimator we need to know the charmed hadron mass, its flight 
vector and the momenta of visible decay products. As stated in ref. |Tl]| the wrong as- 
signment can amount to about 5% for D° as well as for D + and to about 15% for Df. 
In the same paper it was shown that the observed number of DD events is consistent 
with that expected from the charm production cross section, branching fractions, and our 
acceptance. 

The resulting sample of 557 pairs was used for a study of azimuthal angle difference 
A(f) and pseudorapidity gap distributions @]. Since the production characteristics depend 
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only slightly on the nature of the charmed hadrons, we ignore this dependence in the 
present study. 

The simulation of the geometrical acceptance of our apparatus and of the selection 
criteria requires a complex Monte Carlo program. This program generates an uncorre- 
lated pair of charmed particles Pj and Pj. The particles decay into the observed channels, 
and we calculate the acceptance for each combination of decays, each decay being 
generated from phase-space distribution. The distributions of the xf and of the trans- 



verse momenta of the Pi and Pj, as well as the lifetime of A+JT2J are taken from this 



experiment, and the lifetimes of the charmed mesons are taken from RPP[^]. The gener- 
ated pair of charmed particles is mixed with tracks from an event randomly chosen from 
an interaction-trigger sample recorded with our apparatus, thus faking a real event. All 
tracks from such an event are subsequently traced through the apparatus and subjected 
to the same cuts as those made in our analysis. The acceptance vanishes in the back- 
ward centre-of-mass hemisphere, otherwise it amounts to (1 — 5)% depending on the decay 
channel. The acceptance depends only slightly on p?, with a small decrease for increasing 
p T . All results presented in this letter have been corrected for the acceptance. 



3 Momentum estimator 

The momentum estimator, used in order to account for the unseen decay products, 
was first applied to study differential cross sections || and the correlations of charmed 
hadrons |J in the E653 experiment, where the charm decay vertices were observed in 
nuclear emulsion. The NA32 system of the beam hodoscope and the silicon microstrip 
detectors supplemented with two CCDs yielded a comparable precision of measuring the 
primary and the decay vertex, thus allowing us to find the direction of flight of the charmed 
hadron with a high level of accuracy. The visible decay products were identified and their 
momenta were measured with a high level of accuracy by the NA32 large-angle spec- 
trometer. The estimation of the charmed hadron momentum is based on the assumption 
that in its rest frame the invisible system momentum is perpendicular to the charmed 
hadron laboratory momentum. This assumption, the angle between the charmed hadron 
flight vector, and the visible momentum vector, uniquely fix the boost from the charmed 
hadron rest frame to the laboratory frame, thus allowing calculation of the laboratory 
momentum of the charmed hadron. The momentum estimator was obviously not used for 
fully reconstructed decays in our data sample. 

The systematic error of the momentum estimator has been determined using a special 
procedure IfU}] . The charmed hadrons have been generated according to the distribu- 
tions measured in the NA32 experiment H . Then, the hadron undergoes the multibody 
isotropic decay with momenta chosen according to the phase-space distribution^. Next, 
the momentum of the visible decay products and the charmed hadron flight vector are 
found, and the momentum estimator is applied. Then, various kinematical variables for 
the charmed pair are calculated, using the estimated momentum, and finally they are 

4 We belive this is a good approximation in our case since, the majority of our charmed hadrons are 
pseudoscalar mesons, and the main source of error of the momentum estimator arises from the angular 
distribution of the decay products, rather than a particular choice of their momenta. 
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Table 1 



D+/D" decay 
mode 


Api ab /pi ab 
[%} 


^PT,sum 

[GeV 2 /c 2 ] 


A Vdiff 
[1] 


AM eff 
[ GeV/c 2 ] 


^-^(Tr ) 


15.1 


0.086 


0.229 


0.124 




10.0 


0.078 


0.157 


0.099 




18.7 


0.094 


0.276 


0.143 


(K°)7T + 7T + 7l- 


14.2 


0.087 


0.222 


0.123 




15.8 


0.089 


0.244 


0.137 


(K°7r°)7r + 7r- 


19.6 


0.095 


0.292 


0.148 


(K 7r )7r + 7r+7r- 


14.9 


0.090 


0.230 


0.129 


Average 


15.6 


0.089 


0.238 


0.130 



Table 1. The systematic errors of the momentum estimator for various decay modes (with 
unseen decay products in brackets). The numbers are RMS values of the distributions 
of the estimated minus the generated kinematical variable of the charmed pairs. The 
average errors are calculated using weights proportional to the abundance of the decay 
channel in our double tag sample. 



compared to the generated variables. The errors for various decay modes are collected in 
Table 1. The systematic errors of various kinematical variables are determined for seven 
decay modes, which constitute a fair majority of events in our sample. The distribu- 
tions of the estimated minus the generated kinematical variables of the charmed pairs are 
symmetric and centered at zero. An example of such distributions of the differences are 
plotted in Fig. 1 for the D° — > 7^~7r + (7r°7r°) mode with unseen 7r°'s. This decay channel is 
almost twice as abundant in our double charm sample than any other decay mode listed 
in Table 1. In the second column of Table 1 we also quote the estimation error of the 
laboratory momentum of the charmed meson, relative to the generated one. As expected, 
the error increases with the increasing fraction of energy carried by invisible decay prod- 
ucts. Calculating the weighted average we obtain 15.6% accuracy for our estimate of the 
laboratory momentum. It should be noted that this error has little impact on the recon- 
structed kinematical variables of the charmed pairs, especially on the p\ sum for which 
there is very little dependence on the decay channel. We use the RMS of the distribution 
of the difference between the estimated and the generated variable as the measure of the 
systematic error on each kinematical variable in Table 1, checking that the mean value 
of the difference distribution is negligible in each case. The overall errors are calculated 
using the weights based on the numbers of such double charmed events observed in the 



ACCMOR experiment ||10||. We assume these errors to be valid also for the remaining, less 
abundant, decay modes. 



4 Results 

Using the laboratory momenta estimated in the previous section we calculate the char- 
acteristics of charmed pairs in our sample. The acceptance corrected p\ sum distribution 



4 



Table 2 



Experiment 
Beam 


[GeV] 


No of 
pairs 


[(GeV/c)- 2 ] 


[1] 


[(GeV/c 2 )- 1 ] 


E653 p — emul. 
800 GeV/c [§ 


38.7 


35 






0.53+8:15 


WA75 7r~ — emul. 
350 GeV/c 


25.6 


177 


0.50±0.10 


1.00±0.06 a > 


1 -I 7+0.13 
J-- 1 ' -0.17 


NA32 7T - - Cu 
230 GeV/c 
X 2 /NDF 


20.8 


557 


0.50±0.02 
26.1/18 


0.65±0.02 
12.5/9 


1.39±0.06 
8.1/5 



a ) calculated from the measured mean value of the \ydiff\ distribution. 



Table 2. The results of the maximum-likelihood fits to the NA32 data and the corre- 
sponding results of other experiments. 



is shown in Fig. 2. We have fitted the p\ sum distribution with 



do 



,,2 



ne—Pr^ (2) 

u PT,sum 

using the maximum-likelihood method. We have done a similar fit for the distribution of 
the rapidity difference of the charmed pairs 

do - V -W 

~e . (3) 



dy di ff 

In Fig. 3 we plot the absolute value of the rapidity difference \ydiff\- Finally, the distri- 
bution of the effective mass M e ff of charmed particles is fitted to 

<IfT ^ e"«« (M ef f-M ) ^ ( 4 ) 



dM, 



ff 



where M =2m£,. The M e ff distribution is shown in Fig. 4. The fitted parameters, their 
statistical errors, the x 2 values, and the NDF (number of degrees of freedom) are collected 
in Table 2. The results are compared with the results of other charm hadroproduction 
experiments. As can be seen, the yafj an d M e ff distributions become broader with 
increasing y/s, in accordance with QCD predictions. 
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Table 3 



Result 


[GeV 2 /c 2 ] 


< \Vdiff \ > 
[1] 


<M eff > 
[GeV/c 2 ] 




0.434 


0.631 


3.851 


al + h.f. 


0.396 


0.559 


4.467 


NA32 


1.98±0.11±0.09 


0.54±0.02±0.24 


4.45±0.03±0.13 



Table 3. The O(af) QCD predictions (first row) supplemented by the hard charmed quark 
fragmentation (second row) compared to the NA32 experimental results (bottom row). 
The first error is the statistical error of the mean value of the distribution. The second 
error is the systematic error due to the uncertainties of the momentum estimator. 



5 Comparison with O(af) QCD predictions 

The next-to-leading order QCD predictions for the charmed pairs production, calcu- 
lated using the program from ref.||, are collected in Table 3. The mean values of the 
distributions^ were calculated using HMRSB|L6| parametrization of the nucleon and the 
central set SMRS[T7] parametrization of the pion structure functions^. The default value 
of A|^-=122 MeV was used, and the hr and \if scales were chosen to be m c and 2m c , 
respectively, where m c =1.5 GeV/c 2 . The first row in Table 3 represents the results for 
bare charmed quarks with xp>0. As the fragmentation generally leads to the softening 
of the distributions, we use the hard fragmentation function 

f(z) ~ 5(z-l), (5) 

where z = Ed/E c , is the fraction of charmed quark energy transferred to the charmed 
hadron. The results {xp > for both charmed hadrons) are presented in the second row 
of Table 3. Our results, this time determined as the acceptance weighted mean values of 
each distribution, are shown in the last row of Table 3. First errors are the statistical 
errors of the mean values. The systematic errors are taken from the bottom row of Table 
1. 

For the M e ff and Udiff distributions there is a good agreement with the 0(a 3 ) 
QCD supplemented by the hard hadronization of the charmed quark. On the other 



5 The predictions displayed in Table 3 are the mean values calculated using the program from ref.Q, 
and its (signed) weights. These mean values are consistent with the mean values determined from the 
histograms produced by the program, with the maximum number of bins allowed. 

6 It should be noted that the results are not very sensitive to the structure functions, e.g., taking other 
(SMRS) pion structure functions results in ±0.01 GeV 2 /c 2 , ±0.02, and ±0.02 GeV/c 2 differences for the 
mean values of the Pr sum , Udiff and M e ff distributions of the charmed hadron pairs, respectively. 
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hand, the experimental value of p? sum * s much larger than the predicted value. Let us 
recall that this result depends weakly on the uncertainties associated with the momentum 
estimator (cf. Table 1). It should also be mentioned that <p\ sum > = (1.52±0.34) GeV 2 /c 2 
for 20 pairs of fully reconstructed charmed particles in this sample[0. Moreover, the 
WA75 collaboration |15| measured <p 2 T ^ um > = ( 2 - QQ -oil) GeV 2 /c 2 . All these results are 
qualitatively consistent with the relatively flat A0 distributions described previously 

To conclude, the perturbative next-to-leading order QCD calculation, whilst quite 
successful in determining total cross sections, single particle spectra and "longitudinal" 
correlation distributions of charmed pairs (ydiff or M e ff), is not able to reproduce the 
"transverse" correlation distributions (pt sum an d As discussed in refs.|7|, [18|] a fairly 
large intrinsic (or "primordial") transverse momentum of the colliding quarks, comparable 
with the measured mean value of p\ sum , is needed to reproduce the experimental results. 
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Figure captions: 



Fig.l. Plots of the estimated minus the generated kinematic variables of charmed 
pairs like: the Vr^ um (b); the ydiff (c); and the M e ff (d); with one decay vertex being 
D° — > A' _ 7r + (7r 7r ) mode (7r°7r° system unseen). The same difference is plotted for the 
laboratory momentum (a) of a single D° decaying as above. 

Fig. 2. The acceptance corrected distribution of the vector sum of transverse momenta 
squared of the charmed hadron pairs Pt sum- The points are experimental data, the solid 
line is the result of the maximum-likelihood fit, with 6 P =(0.50±0.02) (GeV/c)~ 2 . 

Fig. 3. The acceptance corrected distribution of the charmed hadron pairs rapidity 
difference \ydiff\- The points are experimental data, the solid line is the result of the 
maximum-likelihood fit, with a y =0.65±0.02. 

Fig. 4. The acceptance corrected distribution of the charmed hadron pairs invariant 
mass M e ff. The points are experimental data, the solid line is the result of the maximum- 
likelihood fit, with a M =(1.39±0.06) (GeV/c 2 )- 1 . 
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